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SUMMARY 

Results of correlative and design studies for transition location, laminar 
and turbulent boundary-layer parameters, and wake drag for forward swept and 
aft swept wings are presented in this paper. These studies were performed with 
the use of newly developed, improved integral -type boundary-layer and 
transition-prediction methods. Theoretical predictions were compared with 
flight measurements at subsonic and transonic flow conditions for the variable 
aft swept wing F-14 aircraft for which experimental pressure distributions, 
transition locations, and turbulent boundary-layer velocity profiles have been 
measured. Flight data were available at three spanwise stations for several 
values of sweep, freestream unit Reynolds number, Mach numbers, and lift 
coefficients. Theory /experiment correlations indicate excellent agreement for 
both transition location and turbulent boundary-layer parameters. The results 
of parametric studies carried out during the design of a laminar glove for the 
forward swept wing X-29 aircraft are also presented. These studies included 
the effects of a spanwise pressure gradient on transition location and wake 
drag for several values of freestream Reynolds numbers at a freestream Mach 
number of 0.9. 


INTRODUCTION 

The rise in jet fuel prices during the last decade, and uncertainties of 
the fuel costs for the future have been the cause of continued worldwide 
interest in improving military and commercial aircraft fuel efficiency. 
Additionally, improvements in the range of military aircraft are additional 
factors for achieving aircraft fuel efficiency and drag reductions. The most 
promising aerodynamic means of achieving significant gains in aircraft fuel 
efficiency is to design wings that provide extensive regions of laminar flow 
(refs. 1 and 2). Laminarization of the boundary layer on the surface of an 
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aircraft wing can be accomplished by the use of concepts such as Natural 
Laminar Flow (NLF), Laminar Flow Control (LFC), and a combination of NLF and 
LFC which is referred to as Hybrid Laminar Flow Control (HLFC). 

Extensive regions of laminar flow can be maintained on a wing surface with 
the NLF concept by suitably tailoring the spanwise and chordwise pressure 
gradients to limit the amplification of disturbances that trigger the 
transition of the laminar boundary layer to a turbulent one. Flight tests at 
the Dryden Flight Research Facility on a variable sweep TACT F-lll Fighter 
aircraft (ref. 3) with an NLF wing glove have demonstrated that natural laminar 
flow can be maintained over significant areas at transonic speeds. Recent 
flight tests of a laminar clean-up glove (a thin, smooth fiberglass coating) on 
an F-14 aircraft with a variable-sweep wing (ref. 4), equipped with 
instrumentation for the accurate determination of transition, also revealed 
similar phenomena. Correlative studies were performed during the present study 
for the F-14 aircraft laminar wing, clean-up glove using the methods of 
references 5 and 6 extended to account for spanwise pressure gradient 
effects. Results of these studies are presented in this paper in the form of 
theory /experiment correlations for transition location and turbulent boundary- 
layer parameters. In addition, comparisons of the hot-film experimental data 
on the F-14 aircraft's wing glove with theoretical predictions are described. 

The design of either forward- or aft-swept laminar flow wings at subsonic, 
transonic, or supersonic speeds, using NLF, HLFC, or LFC concepts requires the 
use of reliable theoretical computational methods. A number of finite- 
difference boundary-layer methods, stability methods, and full Navier-Stokes 
equation solvers are available; however, these methods are not well suited for 
the large number of calculations required for parametric wing design and 
optimization studies. Several new theoretical integral -boundary-layer methods 
(references 5, 6, 7, and 8) have small computer execution times and are very 
simple to use. These codes have the capability of predicting the laminar, 
transitional, and separating turbulent boundary layers for applications at 
subsonic, transonic, and supersonic Mach numbers. 

The design studies reported in this paper involve the design of an NLF 
glove for the forward swept X-29 wing. The methods of references 5, 6, 7, and 
8 have been modified for the present theoretical investigations to allow 
consideration of the effects of spanwise pressure gradient. Results are 
presented utilizing these modified codes for various spanwise pressure 
gradients and a prescribed streamwise pressure distribution to determine the 
extent of laminar flow on the corresponding X-29 wing glove. The wing 
sections' upper-surface geometrical changes required to obtain the prescribed 
pressure distribution are also given. Furthermore, the present design- 
optimization studies show qualitatively the procedure for optimizing values of 
spanwise gradients so as to minimize the wake drag with due consideration of 
the effect of imbedded shocks on the boundary layer. Finally a comparison of 
the extent of laminar flow on "equivalent" forward-swept and aft-swept wing 
configurations is made. 
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NOMENCLATURE 


AR 

ASW 

c 


^ D wake 

C P 


sonic 


AC P or AC p 
K 1 c 



e 

FSW 

H 

L, M, N 
K, X 


L.E. 

M 

N.F. 

P 


aspect ratio 
aft swept wing 

chord 2 

C L 

induced drag coefficient - — — 
wake drag coefficient 

pressure coefficient or specific heat of air at constant 
pressure 

pressure coefficient corresponding to local Mach number 
of 1.0 

defined in figure 18 

indicator of either growth or decay of disturbance for U 
component of velocity in transformed plane 

downstream rate of propagation of disturbance for U 
component of velocity in transformed plane 

indicator of either growth or decay of disturbance for W 
component of velocity in transformed plane 

spanwise rate of propagation of disturbance of W component 
of velocity in the transformed plane 

wing efficiency factor for induced drag 

forward swept wing 

* 

X 

form factor = 

groups of dimensionless parameters in transformed planes 
presenting the effects of suction, pressure gradient and 
curvature of velocity profiles. Each symbol defined in 
text. 

leading edge 

Mach number (used only with subscripts) 
transition or instability due to normal flow 
static pressure 
Prandtl number 
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q 


Qx 

or 

Q z 


R 


R T,n 

R inst 

R tran 

s,S,z 


S.F. 


T 

T.E. 

u,v,w 

u,v,w 

x,y 

x, Y 


S' 


1 


f 2 

♦l 


(y) 


dynamic head 

dimensionless parameter representing effects of suction and 
pressure gradient on the shape of velocity profile in 
transformed plane 


n = + VLs . _•?.* 

Q * v 0 dx v 0 • Q z - -57 dr + — 


I 


freestream Reynolds number based on chord = 



value of Ry >n at neutral point 
value of Ry n at transition 

curvilinear coordinates in physical plane; s is normal to 
local sweep lines on wing surface, 5 normal to wing 
surface, and z normal to s-5 plane in spanwise direction 

transition or instability due to spanwise flow 

temperature, °R 

trailing edge 

velocity components in physical planes in s, 5 and z 
directions, respectively 

velocity components in transformed planes in X, Y, and Z 
directions, respectively « 

wing section coordinates in streamwise direction 

Stewartson's transformed coordinates 


s a P 

X = / «/^ds; 



i 


1 


disturbance stream function for transformed plane X-Y 

disturbance stream function for transformed plane Z-Y 

disturbance amplitude function for X direction in 
transformed plane 
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a l 

a 2 


M 

a 


v 

0 


• , X 



0 


P.s 


0 


P.Z 


0 


p.n 


disturbance amplitude function for Z direction in 
transformed plane 

2tt/( wavelength of disturbance in transformed plane in X 
di recti on) 


2n/( wavelength of disturbance in transformed plane in Z- 
di recti on) 

dynamic viscosity 

proportionality constant for viscosity-temperature 
relationship 

kinematic viscosity 


momentum thickness in physical s, £ plane for 2-D or 3-D 

oo ^ 

normal direction 0 = / (1 - xp) d? 

o ”e e e 

momentum thickness in transformed X, Y plane for 2-D or 3-D 


/ 0 - 
o e e 


momentum thickness in transformed Z, Y plane 

= /■/(!-/) dY 

e e 


p,s 


r -ey- ri - 

J D U ' Uj * 


o e e 


/ P ’VS- (» - 5-)* 


o e 


/ p,n ^ 


e e 


<5 greater of 6„ _ or 6 
p,n a p,s p,z 
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6 


boundary-layer thickness In physical plane 


* 

6 


- / 6 a 


pu 


P u 
p e e 


)dC 


boundary-layer thickness in transformed plane 


* 


6 


T,X 


5 t,z 


** 


T>X 


** 


T,Z 


A 

Y 

n 


Subscripts 


o 

e,s 

e,z 

aw 

m 

N 

P 


displacement thickness in transformed X, Y plane for 
2-D or 3-D 

(i - 1- )dv 

o e 

displacement thickness thickness in transformed 1 , Y plane 

- /* ( i- r )<iV 

o e 

/ T,n (1 - £-) ^ dY 
o e 

/ T,n (1- /) dY 
o e 

S T,n 9 reater of 5 T,X or S T,Z 
local sweep angle in degree 

ratio of specific heats for air 

= 2y/b 

stagnation condition 
edge of boundary layer in s-di recti on 
edge of boundary layer in z-direction 
adiabatic wall condition 

mean or average value over boundary-layer thickness 
normal section 
physical plane 
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p,s 


P>Z 


p»n 

s 

T 

T,X 

T,Z 

w 

X 

Z 


L.E. 


T.E. 

6 

0 


00 


s 


physical plane in s-di recti on 
physical plane in z-di recti on 
either p,s or p,z 

direction normal to local sweep lines 

transformed plane 

transformed plane X-di recti on 

transformed plane in Z-direction 

conditions at wall 

transformed X-di recti on 

transformed Z-direction 

conditions at leading edge 

conditions at trailing edge 

parameter based on boundary-layer thickness 

parameter based on momentum thickness 

freestream value 

measured in streamwise direction 


In addition to the above list of symbols and subscripts, symbols and 
subscripts are also defined at appropriate places in the text. 


DESCRIPTION OF THEORY 

A description of the theory for the integral methods used to compute the 
laminar boundary layer, transition, and the turbulent boundary layer on a swept 
wing is given in references 5, 7, and 8. These theoretical methods were 
developed for infinite-span swept wings at subsonic and transonic speeds, and 
for two-dimensional airfoils and bodies of revolution at supersonic speeds. 
However, they did not include the effects of spanwise pressure gradient and the 
resulting effects on the location of transition. The boundary-layer "wash-out" 
phenomena are present on finite aft swept wing configurations with a given 
taper and twist, whereas the boundary-layer "wash-in" phenomena are present on 
configurations with forward sweep. The above theoretical methods of references 
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5, 6, and 7 have been extended in the present paper to account for the effects 
of spanwise pressure gradients on laminar boundary-layer transition. At the 
present time, the effects of suction on the extended method are not considered. 

The theoretical programs, which have been developed and used in the 
present theoretical investigations of the wing glove of the F-14 aircraft and 
for the parametric design studies of the forward swept wing of the X-29 
aircraft with a laminar flow glove, are classified and enumerated as follows: 

1. Potential/Viscous Design and Analysis 

a. Perturbation Method of Characteristics for inverse design at transonic 
Mach numbers. 

b. Bauer-Garabedian-Korn Program for transonic analysis. 

c. NASA-Lockheed Mul ti -component Airfoil Program for subsonic analysis. 

2. Boundary-Layer Analysis 

a. Integral Compressible Laminar Boundary-Layer Method with arbitrary 
chordwise pressure distribution and spanwise pressure gradient, sweep 
and suction at subsonic through high supersonic speeds. 

b. Short bubble and reattachment criteria. 

c. Instability and transition prediction due to leading-edge-normal and 
spanwise flow including the effects of spanwise pressure gradient. 

d. Separating turbulent boundary-layer method from subsonic to transonic 
speeds. 

3. Method for computing C n for an infinitely swept wing from subsonic 

u wake 

through supersonic speeds. 

4. Modified strip method to account for the effects of taper, sweep, and twist 
on wing-section characteristics. 

The theoretical development ~of the laminar boundary-layer method for an 
infinitely swept wing will be briefly described. In addition, the methodology 
for determining the effects of spanwise pressure gradients on laminar boundary- 
layer growth and on transition location will be outlined. 


Governing Equations 

The usual governing equations for a compressible hydrodynamic boundary 
layer on an infinitely swept wing for the coordinate axis system shown in 
figure 1 are the continuity, streamwise momentum in s-direction, and momentum 
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c. ^ 



in the 5 and z directions. The coordinate axis of figure 1 is chosen for the 
purpose of facilitating the application of the present theory to the design of 
the wing sections for an arbitrary wing planform. Moreover, the determination 
of the effect of the spanwise pressure gradient on the extent of laminar flow 
on swept wings by the present methods is greatly facilitated by the use of the 
coordinate axis system selected. Thus, 


+ 3 ^ pV .) = 0 Continuity (la) 

as 85 


P u 




aP . a t au> 
as 85 85 ' 


->• (s-direction 
momentum) 


(lb) 


aP 

n 


o 


(Normal Momentum) 


pu 


aw 

as 


+ 




aw ., 

9? 


(z-di rection 
spanwi se 
momentum) 


(lc) 

(l d) 


The above equations (1) contain the terms consisting of variable physical 
properties, such as density p and dynamic viscosity p. Variations of these 
properties across the boundary layer as well as along the flow direction are 
not negligible for laminar boundary-layer flow at transonic speeds and are 
accounted for. At supersonic speeds, variations in the physical properties are 
very large which means that there is a strong coupling between solutions of the 
hydrodynamic and thermal boundary-layer equations. 


In order to simplify the governing equations (1) for solution by integral 
techniques while maintaining realistic computational results for the 
hydrodynamic and thermal boundary layers, the use of the following Stewartson s 
transformation is made: 


X = 


s a a P c 

f 6 G . 

I a p - ds 

0 a o o 


Y = 


e ? T e 

IT l 


a o p o 0 


( 2 ) 
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where a is the proportionally constant for the viscosity-temperature 
relationship, namely 



T 

T~ 

o 


a = 



T o + 


T aw + 


198.6 

198.6 


(3a) 

(3b) 



(1 + /Pr ^ M^) 

(1 + ^ Hg) 


(3c) 


Furthermore, if it is assumed that the viscous flow on the wing can be 
divided into a finite number of suitably oriented strips on its surface, and 
that each infinitesimal control volume of the strip is a portion of an 
infinite-span swept wing, then the governing equations of motion in the 
physical and transformed planes can be further simplified. The transformed 
boundary-layer equations using Stewartson's transformations, for infinite-span 
swept wing conditions, which can be derived from equations (1) with adiabatic 
wall temperature, are written as. 


3U aV 
9X 3Y 


= 0 


U 


U 


3U w 9U _ 

ax v 3Y 


aw v aw 

ax ¥ aY 


3 f!s 

J e dX 


+ v. 




(4a) 

(4b) 

(4c) 


The relationships between the velocities within the boundary layer and the 
pressure gradients in the physical and transformed planes can be written 
(Pr « 1) as 


a a 

U = — — • u V — — — v ( for adi abati c wal 1 ) ; W * w ( 5a ) 


dU . dM t 5 , 

e _ 1 , e ri , y -1 m2> 1 4 

Wr ~ o a o di“ f 1 + ~T M e^ 


(5b) 
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while the applicable boundary conditions in the transformed planes are given 
by. 


@Y = 0*U = W = 0, 


V = V = - — v 
s aw 
e 


(6a) 


* Y ■ 5 T,X * u ' u e ■ a o M e’ It ■ 0 


(6b) 


and 


@ Y - 6 t>z . W = W e = w e , 


3W 3 W 
aY 


= 0 


aY 


(6c) 


The compatibility conditions at the surface of the wing in the transformed 
planes can be written as, 


au> 


dU. 


" = v,(lv-L ■ H 


•» 2 u 


C (0 U'v 

s^aY^w “ u e dX V. v 2' w 

9 T W 


(7a) 


» Y ■ 0 . - V = v (^) 

s 3Y Z w 0 BY" 3 w 


(7b) 


and in the absence of suction for flows with spanwise pressure gradients, 


1 aP = ^e a^W 
p 9z Y=0 a o V ° aY 2 


I Y=0 


(7c) 


The integral equations in the transformed planes in the leading-edge normal 
direction X and in the spanwise direction Z can be derived by making use of 
equations (4), boundary conditions (6), and wall compatibility conditions 
equations (7), and by use of Leibnitz's rule. 

Momentum Integral Equation in X-Direction, Infinite Swept Ming Conditions 


U e dX ('dlf /dx) = 2 [ L " + ^ ~ s e ] 

C A 


(8) 
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Momentum Integral Equation in Z-Di recti on, Infinite Swept Wing Conditions 


W [“e “e L T ' n l (1 - 

0 e e 


-aw 


v oUy)y= 0 ‘ V s W e (9) 


where 6y n is greater of 5-j- z or S-j- x - 


The symbols L, K, and S Q appearing in the above equation (8) are defined 
in equation (20) for transformed X and Z coordinates. 

The velocity profile assumptions for the flow in the X- and Z-directions 
in the transformed planes are given by the following equations, 


where 


IT = A l n x + Vx + A 3 n X + Vx 
e 



(10a) 


and 


- 2 + Q x / 6 ; ~ - 0.5Q X J 

A 3 = -2 - A 2 ; A 4 = 3 - Aj 


(10b) 






(ID 


where, n 7 = t , 

Z 6 T,Z 


The numerical solutions of equations (8) and (9) are performed by the 
Euler method with repeated iterations. These solutions give the shape of 
velocity profiles jj- and jjjj- in the transformed planes under the assumption of 

e e 

infinite-span swept-wing conditions. 
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EFFECTS OF BOUNDARY-LAYER HASH -IN/WASH -OUT DUE TO SPANWISE 
PRESSURE GRADIENTS ON TRANSITION ON FINITE, SWEPT, 
TAPERED, AND TWISTED WING 


The results of theory /experiment correlations for several wings using both 
wind tunnel and flight-test data have indicated that for sweeps greater than 
approximately 25°, and freestream Reynolds numbers greater than 20 million 
transition is triggered by spanwise flow in the z direction. Thus, for the 
purpose of designing laminar-flow wings and gloves at transonic and supersonic 
speeds at Reynolds numbers and sweeps of practical interest, it is necessary to 
consider the effects of boundary-layer wash-in/wash-out due to spanwise 
pressure gradients. This is especially important for wings of supersonic 
aircraft where the sweep of the leading edge is usually in excess of 45°. In 
order to account for the effects of boundary-layer wash-in/wash-out on finite 
swept wings the approach shown schematically in figure 2 is used. This figure 
shows the infinitesimal control volume composed of surfaces of a trapezoid. 

The coordinate system is the orthogonal s, e; and z axes where the s 
coordinate is normal to the local sweep line, z-coordinate is normal to s in 
the plane of wing, and the ^-coordinate is normal to the s-z plane. 
and are constant s/c lines where c is the local chord in s-direction 

and and are parallel to s-direction. 


The momentum theorem is applied in the z-di recti on, which states that the 
algebraic sum of the forces acting on the control volume ABCD - A^B^C^D-^ in 
the z-di recti on is equal to the rate of change of momentum in control volume 
plus the net flux of momentum in the z-di recti on across the control surfaces of 
the control volume. For steady mean flow conditions an equation can be derived 
for the boundary-! ayer wash-in/wash-out due to spanwise pressure gradient 


effects on velocity profile — in the z-di recti on. This equation in the 

w 

e 

physical coordinate system for transformed boundary-layer quantities can be 
written as: o 


dC. 


k (0 T,z } N = { az [{Q l,l ) o ~2 (1 + °* 2 M e,z )0 ‘ 5] 

00 

+ W- IT (1 * °' 2 M e,s>°' 5 [ i W ljL (5 T,X >0 (0 - 7 + °- 00833 V 


OO 00 


e,s 


dA. 


H Y 

+ < 0 - 7 + °- 00833 V Ss (6 T,X } o + 0-00833(6 T>x ) o ^ ] 


w w W e 1 1 n (1 + 0.2 M e,s 2 r 3 , 

— * v n* rr * JT (e T z ) q * °* 235 " — ?.) 


(1 + 0.171 M 

c j S 


(1 + 0.2 M 


2 ) 
e,z 


-3 


(1 + 0.2 M 2 ) 2,5 / (8.547 + 0.711 M 2 ) 

00 6 , Z 


( 12 ) 
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where 


( 0j z )u = Transformed momentum thickness in Z-direction 

’ which is modified by the boundary-layer wash -in/ 

wash-out phenomena due to spanwise pressure 
gradients on finite-swept and tapered wing. 


<V»o ■ 


Transformed momentum thickness in Z-direction for 
infinitely swept wing assumption conditions. 


U T J X ) o J 0 



Dimensionless pressure gradient parameter in 
X-di recti on in transformed plane for infinite-span, 
swept-wing assumption conditions. 


Transformed boundary -layer thickness in X-direction 
for infinite-span swept-wing assumption condition. 


Tp + 198.6 <j ue tQ Sutherland 

T + 198.6 theory of viscosity, 
aw 

approach is used for the purpose of calculating (6y ^ and 
the velocity profiles in the transformed plane in the spanwise direction in 
order to account for the boundary-layer wash-in or wash-out phenomena. 

Define the quantities. 


and. 


where: 

( 6f ^ ) |g = Transformed boundary-layer thickness in Z direction which is 

’ modified by the boundary-layer wash-in/wash-out phenoena due 

to spanwise pressure gradients on finite-swept and tapered 
wing. 

If a fourth order velocity profile is assumed in the spanwise Z direction 
in the transformed plane, then the relation between A, and a can be derived 

Zi. 2 

as the following polynominal: N 


2 1 dW e 

V? = ( 6 y 7 )y ( — ~ ) 

Z N T ’ Z N v o d 2 


2 1 dW e 

e 7 = (0 y 7 ) w (— *£) 

Z N T,z N v 0 dZ 


(13) 


( 6y y) n = 

T,X 0 


,1.5 


w / aw 


The following 


180 



where. 


h = P 1 + P ? S 7 + P ^ 2 7 + P A&7 + P q4 + P fi*7 

Z N 1 2 Z N 3 Z N 4 Z N 5 Z N 6 Z N 


P x = 0.032015 
P 2 = 73.156 
P 3 = 10.587 
P 4 = -129.6 
P 5 = 4669.6 
P 6 = 43865.0 

thus, (« T>Z ) N can be written as, U T>Z ) N = (0 T Z ) N / (e z /A z ) 1/2 . 


w 

The velocity profile y- in the transformed plane can therefore be calculated 

e 

from the following assumed fourth order polynomial given by: 


_W _ R n . 

W B l n Z 
e 


B^n 


2" Z 


3 
3"Z 


B^n 


+ B 


4 Z 


where. 


n, - 


(fi T,Z 

Bo 


Z N 1 

and Bi - 2 + — y — ; B „ = - A 

c Z N 

-2 - B 2 ; B 4 = 3 - B x 


THERMAL BOUNDARY-LAYER SOLUTION 

Temperature variations within the boundary layer at subsonic speeds are 
small, hence the effects of temperature variations within the boundary layer on 
instability and transition at subsonic speeds are neglected. However, at high 
transonic and supersonic speeds, the heat generated by friction and adiabatic 
compression is quite significant. These phenomena give rise to large 
temperature variations across the boundary layer, and it has been found from 
computational experiments that these large temperature variations affect the 
following equations for computing temperature profiles; they can be derived 
using the procedures described earlier. 
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(14a) 


w 


'Y=6 


T ,X 


w 


aw 


e Y = 6 


T,X 


(tt) - ^ 6 


2gC T 

3 p 0 


t )2 


T - T. 


w 


"Y=6 


T,Z 


w 


aw 


Y=6 


T,Z 


V 


- /Pr 


2gCT 
3 p o 


e 


(14b) 


For determining the temperature profile at a given location on the wing 
use is made of either equations (14a) or (14b) depending upon whether is 


in the 


U W 

larger than <5 T 7 or vice versa. The velocity profiles tt- or 

1 * L e e 

transformed planes in X or Z direction, which are required in equations (14a) or 
(14b), are dependent on local pressure gradients in the leading edge normal -s 
or spanwise-z directions due to use of equations (9) or (11), respectively. 

Thus, the effect of pressure gradient and suction is implicitly accounted for 
in the determination of temperature profiles by the present method. 


Criteria for Determining Locations of Neutral Instabilities 
for Velocity Profiles in s and z Directions 

The Orr-Sommerfeld equations in transformed planes for X and Z 
directions can be written as follows: For the transformed coordinate in the X- 

di recti on, 

(U - C 1 )(^' - a^l) “ U YY^1 = oTF " _2tx i + “i*].] 

1 9 T,X 

and for the transformed coordinate in the Z-di recti on, 

(W - C 2 ) ( 4>2 ‘ - <*2^2 ) " w yy* 2 = — [+2 " 2< *i+2 ' + (16) 

2 0 T,X 

The disturbance stream functions for the transformed X and Z directions 
are assumed, as for the physical plane, as the following types: 


ia, (X-C,t) 

*1 = *l(Y) e 



+ X- direction 

(17a) 

ia 2 (Z-C 2 t) 



*2 - *2(Y) e 

+ Z-direction 

(17b) 
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where. 


4 >., Y x and 4 > 2 ( y ) = disturbance amplitude function 

1 [ ’ for X and Z direction, respectively 

in transformed planes 

a, and a ? = 2ir/(disturbance wavelength in transformed 

planes in X- and Z-di recti on, respectively) 


X = transformed coordinate for s-di recti on 
Z = transformed coordinate for z-direction 


C 


1 




C 


2 




complex in general 


C, or Cp = 0 
> 0 


< 0 


corresponds to neutral instability 
corresponds to amplification of disturbance 
corresponds to decay of disturbance 


The generalized solution between the dimensionless pressure gradient 
parameters and "equivalent" Reynolds number based on momentum thickness for the 
neutral stability of laminar boundary layers in either s or z direction is 
shown in figure 3. The curve shown in this figure is derived from solutions of 
the Orr-Sommerfeld equations, equations (15) and (16), in the transformed X-Y 
plane or Z-Y plane for appropriate fourth-order polynomial velocity profiles. 
The effects of pressure gradient in the s-di recti on or z-direction, local Mach 
number, temperature profile, and suction have been accounted for in the 
derivation of this curve. The abscissa of figure 3 is given as the zero 
suction ( S = 0) version of 



1 


( 1 + 0.2 


13m! 


[i - 


(M + S 


M 


0 


e,z 


TT ^ 1 ne, 

or 
L, 
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where 


m 


= f / 


n 0 


vd? 


greater of „ or 6 , 

3 P,s p,z 




/ 2 ' 
M 

e,s 

M 

e,z 


1.5 T 

* T 


r^ s 

e,z 


( 19 ) 


In equation (19), the term (y ) or (y ) for temperature profile are 

e,s e,z 

substituted from equations (14a) or (14b). 

It should be noted that because of the choice of the coordinate axis 
systems shown in figures 1 and 2, it is possible to use the same criteria shown 
in figure 3 for determining the neutral stability point for disturbances in 
both the normal and spanwise directions. 


Boundary -Layer Transition Criteria 


The criteria used in the present method for determining the transition 
loc atio n due to the amplification of disturbances either in the normal or 
spanwise directions are shown in figure 4. The parametric curves of figure 4 
are valid for arbitrary specified pressure gradients - both in the normal as 
well as spanwise directions. The independent and dependent variables shown in 
this figure contain groups of several hydrodynamic and thermal boundary-layer 
parameters which are significant during the transition process. In deriving 
the curve for low values of freestream turbulence (less than .05 percent), the 
use of experimental data obtained in flight or low-turbulence wind tunnels is 
made and this information is put in the transformed planes. The dimensionless 
parameters shown in figure 4 have been derived from dimensional analyses and by 
making use of the laws of dynamic similarity. These physical dimensionless 
parameters consist of the following: 
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In deriving the curves for several higher values of freestream turbulence 
greater than 0.05% shown in figure 4, use has been made of information 
presented in figure 5 for a zero pressure gradient and for values of Mach 
number from low subsonic to high supersonic (up to = 4.0). It has been 
hypothesized in constructing the curves of figure 4 for values of freestream 
turbulence greater than 0.05 percent, that the effect of pressure gradient for 
freestream turbulence greater than 0.05 is similar to low values of freestream 
turbulence of 0.05. Curves for freestream turbulence higher than 0.05 percent, 
shown in figure 4, have been constructed by making use of information in figure 
5 and using the above hypothesis. 

It needs to be pointed out here that the curves of figure 5 can also be 
used for other types of disturbances, e.g., noise generated due to the 
turbulent boundary layer on the walls of the test section of supersonic wind 
tunnels. This objective can be accomplished after the application of 
appropriate calibration procedures. 


Turbulent Boundary Layer 

It has been known from experimental data on sections of NLF wings that if 
turbulent separation is present upstream of x/c k 0.95, then the design 
pressure distribution necessary to maintain laminar flow is degraded. 
Additionally, the wake drag will be quite high due to presence of the separated 
flow region. In the present paper turbulent boundary-layer separation is 
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accounted for by the method developed in reference 8. This turbulent boundary- 
layer method has already been successfully applied in the design of NLF and 
HLFC wing sections and for the prediction of wake drag. 

To compute the turbulent boundary layer from the location of transition at 
various spanwise stations to the wing trailing edge, the method of reference 8 
was modified to account for effect of spanwise flow assuming infinite swept 
wing conditions. The important objectives for the computations of the 
turbulent boundary layer are: (1) to see if turbulent separation is present 

downstream of the transition location, and (2) to determine the turbulent 
boundary-layer quantities at the wing trailing edge. 

The theoretical method of reference 8 has been developed by analyzing the 
mean and fluctuating components of velocity profile data obtained by the use of 
a laser velocimeter for several airfoils with trailing-edge seperation. These 
experimental data indicated that as the turbulent boundary layer approached the 
separation location, the fluctuating velocity components become of the same 
order of magnitude as the mean velocity components within the turbulent 
boundary layer. This observation suggested the retention of the fluctuating- 
velocity components in the governing equations of motions as was done in 
reference 8. This method has been found to be extremely reliable in predicting 
turbulent boundary-layer quantities, including separation, on wing sections 
with thickness ratios in the range of 0.09 < t/c < 0.22 and Reynolds numbers 
from 2 to 30 million. 


Computations of Cn for Ming Sections 
u wake 

The derivation of the expression for Cq for an infinitely swept wing 

wake 

is given in Appendix A. This expression, which is valid from subsonic to 
supersonic speeds, is given as 
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where the exponent p is given by 

p = [ h t.E . + 5 + W-l)M^]/2.0 (21) 

The C °wake f ° r a fin1te ta P ered s we Pt wing can be computed from equation 
(21) and by making use of modified strip theory. 
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APPLICATIONS OF THE PRESENT THEORETICAL 1CTHODS 


The theoretical methods described in previous sections are utilized to 
perform two tasks in the following section, namely. 


(i) 


The predition of transition location and turbulent boundary-1 ayer 
parameters on the "clean-up" glove of the F-14 aircraft and 
correlation with experimental data. 


(ii) Performing parametric studies for the design of a laminar glove for 
the forward swept wing of X-29 aircraft. 


Correlation with Experimental Data on the Glove of F-14 Aircraft 

The planform and instrumentation for the gloved F-14 wing with 20 sweep 
is shown in figure 6. With the "clean-up" glove installed the sweep of the 
wing leading edge was varied from A L>E> * 20 to A L>E> - 35 , and the p 

slats were locked in the up position. As shown in figure 6, the three rows of 
surface pressure orifices are located at spanwise stations n = 0.515, O.b/, 

0 83 respectively, and were oriented in a direction parallel to freestream for 
the wing leading edge sweep A, F - 20°. In order to avoid interference 
between hot-film gauges, they L ilre placed at an angle of 30° to rows of surface 
pressure orifices as shown in this figure. The boundary-layer rakes were 
installed at two spanwise locations at x/c = 0.55 in order to measure te 
turbulent boundary-layer velocity profiles. Transition and turbulent boundary- 
layer data were obtained for freestream Mach numbers from 0.6 to 0.8 wing 
leading edge sweeps from 20° to 35°, and a Reynolds number range of 8x10 to 
25xl(j Pressure distributions for the aforementioned three spanwise 

locations were also measured for these Mach numbers and sweeps. These Pressure 
distributions constituted the boundary conditions for computing both the laminar 
and turbulent boundary layers and the transition location on the F-14 g ove. 

Figure 7(a) shows the pressure distributions on the upper surface of the 
qlove at three spanwise locations and hot film traces at station 2 
(n = 0.67) for the wing leading edge sweep of 20 , M m = 0.7 and a freestream 
Reynolds number based on chord of 13 million. In addition to using t e 
pressure distributions as boundary conditions, they are used to evaluate t e 
effect of spanwise pressure gradient on the development of the laminar 
boundary layer and transition location at station 2. Since the experimenta 
pressure data only extended back to x/c = 0.6 theoretically computed pressures 
distributions were used from x/c = 0.6 to trailing edge in order to compute 
the growth of turbulent boundary layers from the transition location to the 
trailing edge. Figure 7(b) shows the computed results for chordwise 
transition and instability location due to normal flow (N.F.) and spanwise 
flow (S.F . ) phenomena. The results are plotted in this figure as x/c of 
instability and transition versus freestream Reynolds number. These results 
indicate that at a Reynolds number of 13 million, the computed transition 
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location is due to N.F. phenomena and occurs at x/c = 0.33. The hot film 
trace shown in figure 7(a) shows that the experimental transition location at 
station 2 is between 0.3 < x/c < 0. 4. Thus, good agreement is obtained 
between the theoretically computed transition location and the hot-film- 
measured transition location. Figure 7(c) shows the results of computations 
for displacement thickness, 6*/c, momentum thickness, e/c, and form factor H, 
plotted against x/c for station 2. Freestream Mach number, Reynolds number, 
and leading edge sweep for the results shown in figure 7(c) are the same as 
for figures 7(a) and 7(b). Experimental data for <s*/c, q/ c , and H at x/c = 
0.55 for station 2 in the turbulent boundary-layer region are also plotted in 
figure 7(c). These results indicate that present theoretical computations for 
the turbulent boundary-layer parameters agree quite well with the 
experi mental ly measured quantities. Furthermore, these results indicate that 
the location of transition predicted by the present theory must be in close 
agreement with experimentally determined transition. It is also interesting 
to notice from figure 7(c) that the slope of e/c versus x/c in the laminar 
flow region is smaller than that in the turbulent region. This large slope of 
e/c versus x/c in the turbulent boundary-layer region will lead, of course, to 
large values of e/c at the wing trailing edge. The wake drag coefficient/ 
C °wake’ iS dlrectly P r °P ort1 ’onal to the momentum thickness, e/c, at the 

trailing edge. 

Figures 8(a) and 8(b) show plots of experimental pressure distributions 
and hot film traces for = 0.801, a leading edge sweep of 20° and freestream 
Reynolds numbers of 11 million and 23 million, respectively. As can be seen 
from the comparison of pressure distributions in figures 8(a) and 8(b), there 
is not any noticeable difference in pressure distributions between Reynolds 
numbers of 11 and 23 million, and hence, the same pressure distributions at 
three spanwise stations were used for computations of transition and 
instability results for the range of Reynolds number shown in figure 8(c). 

These pressure distributions indicate that the shock strength is high, and for 
this reason, there seems to be some discrepancy in computed and measured 
turbulent boundary-layer parameters at x/c = 0.55, as shown in figure 8(e), 
even though the transition location is predicted quite accurately. The present 
shock boundary-layer model is being refined in order to provide better initial 
conditions for turbulent boundary-layer computations downstream of a strong 
shock . 

As seen from figure 8(c), the computed transition location for a 
freestream Reynolds number of 11 million is at x/c ~ 0.5 and caused by N.F. 
phenomena, however, at a freestream Reynolds number of 23 million the computed 
results indicate that transition is caused by S.F. phenomena and occurs at 
x/c = 0.41. Experimental data from the hot” film sensor shown in figure 8(a) 
indicate that at a Reynolds number of 11 million, transition occurs at 
x/c * 0.5. Hot film traces at a Reynolds number of 23 million, which are 
shown in figure 8(b), indicate that transition is occurring, for this 
freestream condition, at x/c ** 0.4. Thus, as seen from comparisons of 
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computed results, shown in figure 8(c), with hot film traces of figures 8(a) 
and 8(b), the present theory reliably predicts both the location of transition 
as well as its mode, i.e., due to N.F. or S.F. phenomena. Figure 8(d) shows a 
plot of computed C n as a function of freestream Reynolds number. As seen 

from this figure, C n decreases as the freestream Reynolds number increases 

u wake 

from 8 to 20 million and above the Reynolds number of 20 million C n 

u wake 

increase with Reynolds number. The Reynolds number of 20 million, where the 

computed C n becomes minimum, corresponds to the initial triggering of 
wake 

transition due to spanwise flow. This is an important phenomena and one that 
can be used to advantage as seen in the design of a proposed laminar glove for 
X-29 aircraft discussed subsequently in this paper. Finally, figure 8(e) 
shows the plots of computed values of 5*/c, e/c and H versus x/c for a 
Reynolds number of 12.4x10^. 

Figure 9(a) shows the experimental pressure distributions and hot film 
traces for a wing-glove leading-edge sweep of 25°, a freestream Mach number of 
0.706 and a Reynolds number of 11 million. By using the experimental pressure 
distributions as boundary conditions (see figure 9(a)), the transition and 
instability locations are computed and the results are shown plotted in figure 
9(b) as a function of freestream Reynolds number. Experimental data for 
transition obtained from hot film traces, shown in figure 9(a), indicate that 
transition for a wing glove leading edge sweep of 25.3°, Reynolds number of 11 
million, and M = 0.706 takes place for 0.4 < x/c < 0.5. Theoretical results 
for these same conditions shown in figure 9(b) also show that transition 
should occur at an x/c between 0.4 to 0.5. 

Now, it is generally believed that if the freestream Mach number is kept 
constant, that increasing wing sweep reduces the extent of laminar flow. 
However, experimental data and calculations for the wing-glove leading-edge 
sweep of 20° and M^ = 0.7, (figures 7(a) and 7(b)) when compared to results 
for a sweep of 25°7 figures 9(a) and 9(b), show that transition is located 
closer to the leading edge for the lower sweep. Furthermore, comparison of 
the calculated results of figure 9(b) with those of figure 7(b) indicates the 
same phenomena. This is an unusual result and shows the ability of the 
present methods to optimize wing section geometries and planforms for 
extensive laminar flow regions at specified freestream conditions. 

Figure 9(c) shows the results of theoretical computations for the laminar 
and turbulent boundary-layer parameters <$*/c, e/c and H plotted versus x/c for 
a leading-edge sweep of 25°, M of 0.704 and a Reynolds number of 
approximately 11 million. For this case, strong shocks are not present at any 
spanwise station, as can be seen from figure 9(a). The computed boundary- 
layer parameters agree quite well in the turbulent boundary-layer region at 
x/c = 0.55. This agreement of computations with experiment in the turbulent 
boundary-layer region is an indication that the computed transition location 
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is fairly close to the experimental transition location. Moreover, these 
results indicate the reliability of the present theoretical turbulent 
boundary-layer method. 

Figure 10(a) shows the experimental pressure distributions and hot film 
traces for a leading-edge sweep of 29.7°, M = 0.704, and a Reynolds number of 
9.5 million. The corresponding experimental data for a Reynolds number of 
19.2 million and leading edge sweep of 31.4° are shown in figure 10(b). 

Figure 10(c) shows theoretical results for instability and transition 
locations for a leading-edge sweep of 30° and M^ = 0.704. These results show 
that for a Reynolds number of 9 million transition takes place at x/c = 0.34 
and is caused by S.F. instability phenomena. Experimental data obtained by 
hot-film sensors, shown in figure 10(a), for a Reynolds number of 9 million 
indicate similar results, namely transition takes place at an x/c between 0.3 
to 0.4. At a Reynolds number of 19.2 million (figure 10(c)), transition due 
to S.F. takes place at x/c ~ 0.175; experimental hot film trace data for the 
same freestream conditions, shown in figure 10(b), also indicate that 
transition occurs between x/c of 0.1 to 0.2. 


Design and Parametric Studies for Laminar Glove of X-29 Aircraft 

It has been observed from theoretical studies and experimental data that 
wing sweep have a pronounced effect on transition location when the freestream 
Mach number and Reynolds number are held constant. It is also believed by 
many, from a conceptual view point, that it is easier to achieve large 
laminar-boundary-layer-flow regions on a forward swept wing (FSW) than on an 
aft swept wing (ASW). This is so because for FSW the local values of sweep 
are lower near the leading edge and the local sweep angle increases from the 
leading edge to the trailing edge; whereas, for ASW, sweep angle is the 
largest near the leading edge and decreases progressively toward the trailing 
edge. Figure 11(a) shows schematically a comparison of planforms for 
"equivalent" FSW and ASW. The word "equivalent" implies that for both FSW and 
ASW, wing area, aspect ratio, taper ratio, and shock locations are 
identical. Figure 11(b) shows the plots of local sweep angle versus x/c for 
FSW and ASW of aspect ratio of 4 and taper ratio of 0.4. The planform and 
variation of sweep angle with x/c shown in figures 11(a) and 11(b) are a rough 
approximation of the X-29 aircraft wing. The point C in figure 11(b), which 
is at the intersection of local sweep versus x/c curves for the "equivalent" 
FSW and ASW, indicates that the local sweep angles for these "equivalent" 
wings become equal at x/c ~ 0.5. Thus the possibilities of maintaining 
laminar flow for x/c < 0.5 on an NLF glove for the FSW of X-29 aircraft are 
much greater than they would be for an "equivalent" ASW with a laminar glove. 


It is possible to realize an additional advantage from forward swept 
wings and this stems from the lower shock-wave drag. If the NLF design 
pressure distribution is such that the minimum C D on the upper surface occurs 
aft of x/c ~ 0.5, then the local Mach number at the shock location for FSW 
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will be lower than that ahead of the shock on an "equivalent" ASW for the same 
value of the freestream Mach number. This is so because the local sweep angle 
for FSW is higher than that for an equivalent ASW downstream of point C 
(x/c ~ 0.5) as can be seen from figure 11(b). These phenomena would result in 
lower shock wave drag for FSW for the NLF design glove of X-29 aircraft, and 
moreover, the susceptibility of the NLF design for FSW to turbulent separation 
is also diminished. This is because the initial turbulent boundary -layer 
thickness at the shock location would be lower for the FSW than for the 
equivalent ASW. 

The above discussion suggests that better prospects exist for overall 
lower total drag for an NLF concept with FSW than one with an equivalent 
ASW. The superior drag peformance due to application of NLF (or HLFC) 
concepts to a FSW arise due to the possibility of achieving larger chordwise 
runs of laminar boundary -1 ayer flow, lower shock wave drag, and less 
susceptibility to turbulent separation drag. 

Figure 12 shows the planform of the X-29 aircraft. The canard and wing ^ 
are in the same plane. The wing has a leading-edge sweep of approximately 29 
and the trai ling-edge sweep is 45°. The chord of the first flap extends from 
x/c = 0.75 to x/c = 0.9 and that of the second flap extends from x/c = 0.9 to 
x/c = 1.0. The wing leading-edge break is at n = 0.25, and the tip of the 
canard is at n = 0.49. Experimental pressure distributions on the canard 
suggest that boundary -layer flow separation exists on the canard and the 
possibility exist that the separated wake of the canard may prohibit the 
achievement of laminar flow on the wing between the wing root and the midspan 
location. The forward swept wing has associated with it the boundary-layer 
wash-in phenomenon, and for this reason it is believed that a laminar glove 
can be designed so as to obtain extensive laminar flow on the outboard wing 
from n = 0.55 to n = 0.95. 


Simple-sweep methods or modified simple-sweep methods are frequently used 
in the design of airfoil sections for aft swept wings. Different forms of 
modified simple-sweep methods are employed but they give essentially the same 
results. The present modified simple sweep approach for ASW is illustrated in 
figure 13. The geometry of the normal section is arrived at from the geometry 
of the streamwise section by the following expression: 


where. 


_ ^ Z//c ^streamwise section 

( z / c ' 1 Normal section Cos 

A l = local value of sweep 
= a le + (a je - a le ) • (x/c) 
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The pressure distribution is computed for the normal section by airfoil 
programs, such as, Bauer, Garabedian, and Korn program (ref. 9) or Lockheed- 
NASA multi -Component airfoil program (ref. 10). The computations for the 
normal section Mach number and lift coefficient are made using 

°» s effective 

where, 

M = effective component of freesteram Mach number for normal section 
N 

M = streamwise freestream Mach number 
00 

S 

and, A e ff ec - t j ve = A 0 x/c = 0.25 for subsonic and transonic freestream 
Mach numbers 


and 

C. _ C, /Cos^A ff .. 

L^j = effective 

where, 

Ci = lift coefficient for normal section 

N — __ _ 

C|_ = lift coefficient for streamwise section 


The thickness of the normal section for the aft wing section is thus 
greater than that of the streamwise section by the factor of 1 /Cosa. . The 
pressure distribution computed for this normal section (with L 

M = M CosA ff and C, = C. / Cos^A is related to the pressure 
N S l N l S ert 

distribution for the streamwise section by, 

Cp = Cp • Cos A. 

Streamwise Normal L 

The application of the above-mentioned modified sweep approach to the 
computation of the pressure distribution for the forward swept wing of the X- 
29 at n = 0.7, M^ = 0.91 and C^ = 0.351 is shown in figure 14. It is 
apparent from this figure that the theoretical pressure distribution computed 
by the above-mentioned modified simple sweep method for aft swept wings does 
not agree with the experimental data for the FSW of X-29 aircraft. 
Consequently a new modified simple sweep approach was developed for the 
computation of pressure distributions for forward swept wings. In this 
approach, the ordinates of the upper surface of the normal section are kept 
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the same as for the streamwise section but the thickness distribution for the 
normal section is increased by the cosine of the local sweep angle and is 
given by 


(t/c) Local = (t/c) Local / CosA. 

Normal Streamwise L 

Streamwise Section 

The remaining steps for the computation of the streamwise pressure 
distribution are the same as for the aft swept wing. Both procedures are 
illustrated in figure 13. The results of computations of the streamwise 
pressure distributions, using the new modified simple sweep method, for the X- 

29 aircraft for n = 0.7, = 0.91 and = 0.351 are shown in figure 14. As 

seen in this figure, the agreement between the theoretical upper surface 
pressure distribution computed by the new modified simple sweep method and 
experimental data is reasonably good. Hence, the parametric studies for the 
design of the outboard laminar-glove sections are conducted using this 
approach. 

In order to enlarge the extent of laminar flow for the X-29 wing, by the 
installation of a wing glove, it is first necessary to determine the extent of 
laminar flow on the existing X-29 wing. Theoretical computations for 
instability and transition locations on the existing X-29 wing at n = 0.7 are 
plotted in figure 15(a) as a function of freestream Reynolds number. 
Experimental pressure distributions are used for boundary conditions in this 
calculation. The computed wake drag coefficient for the existing wing section 
at n = 0.7 is plotted as a function of freestream Reynolds number in figure 
15(b). The results shown in figures 15(a) and 15(b) illustrate that the 
extent of laminar flow on the existing X-29 wing is quite limited and, 
consequently, the wake drag is high for Reynolds numbers in the range of 8 to 

30 mi 11 i on . 

A number of pressure distributions and corresponding glove geometries 
have been theoretically investigated in the present studies, under the 
assumption of infinite-span swept-wing conditions. It is necessary that the 
glove geometry be compatible with original X-29 wing such that the glove wraps 
the wing surface with continuous first and second derivatives, and that there 
be a smooth merging of the glove with the original wing surface. From these 
considerations, the glove with the geometrical shape designated as XTM2P3K, 
and shown in figure 16 for spanwise station o = 0.7 was developed. The 
geometry of the baseline X-29 for the spanwise station of n = 0.7 is also 
shown in figure 16 for comparison purposes. Figure 17(a) shows the computed 
streamwise pressure distribution for this glove at n = 0.7 for M = 0.905, 

C. = 0.435 and a leading-edge sweep of 30°- The computed pressure 
S 

distribution for = 0.351 for the baseline section is also shown In figure 

17(a) for comparison. It is important to note the shape of C p in figure 

K sonic 
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17(a). For the forward swept tapered wing, the C p . line has an upward 

som c 

slope, whereas for the "equivalent" ASW the C p . line will have a downward 

r somc 

slope. These phenomena suggest, for the specified freestream Mach number and 
value of C p just upstream of the shock, that the FSW will have a smaller value 
of shock wave drag than the ASW. 

The streamwise pressure distribution computed for the glove with a L.E. 
sweep of 25°, and = 0.43 is shown in figure 17(b); since the sweep is 

lower, the design freestream Mach number was lowered to a value of 0.865 (as 
compared to an M^ of 0.905 for sweep of 30°). The reason that the Mach number 
was lowered was to keep the leading-edge normal Mach numbers constant. The 
theoretical transition characteristics for the X-29 glove for leading-edge 
sweeps of 30° and 25° and zero spanwise pressure gradients are shown in figure 
17(c). As seen from figure 17(c), the S.F. instability triggers the 
transition at a Reynolds number of 12 million for a leading edge sweep of 30°, 
and at a Reynolds number of 15 million for a leading edge sweep of 25°. 

Figure 17(d) shows a comparison of the wake drag characteristics for the 
X-29 glove for leading edge sweeps of 30° and 25°. Also plotted in figure 
17(d) is the wake drag of the baseline section for a leading-edge sweep of 
30°. These comparisons indicate the following phenomena, first, the wake drag 
of the glove is lower than the basic X-29 wing section by as much as 17 to 18 
counts at a Reynolds number of 10.5 million, and second, when the leading-edge 
sweep is reduced from 30° to 25°, the freestream Reynolds number at which 

C n becomes minimum is increased from 10.5 million to 16 million. 

D wake _ 

However, it may not be possible to reduce the wing glove leading-edge sweep to 
25° because of stability and control considerations. 

As stated previously in the theoretical section, the methods of 
references 5, 6, and 7 have been extended in the present studies to take into 
account the effects of boundary-layer wash-in/wash-out and spanwise pressure 
gradients on the development of the laminar boundary layer as well as 
instability and transition locations. Thus, it was possible to con duct 
parametric studies on the effects of boundary-layer wash-in phenomena and 
spanwise pressure gradients on the triggering of transition due to S.F. 
phenomena. Proper use of spanwise gradients can delay transition to higher 
Reynolds numbers and result in a wake drag reduction at the cruise Reynolds 
number. Figures 18(a) and 18(b) show the two types of spanwise pressure 
gradients which have been utilized to study the effect of spanwise pressure 
gradients on transition. In figure 18(a), constant spanwise pressure 
gradients are assumed on the upper surface from 0 < x/c <0.6; computations 
have been performed for spanwise gradients, ACp , of 0.05, 0.1 and 0.15. For 

c 

the linearly varying spanwise pressure gradients depicted in figure 18(b), the 
spanwise pressure gradient is a maximum near the leading edge and it decreases 
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linearly to zero at x/c = 0.6. Computational results are presented for 
ACp^ values of 0.05, 0.1, and 0.15. Instabilities and transition 

characteristics due to N.F. and S.F. and wake drag characteristics for the 
various spanwise pressure gradients shown in figures 18(a) and 18(b) have been 
determined. A few representative results are shown in figures 19 and 20. 
Figure 19 shows the plots of N.F. and S.F. transition characteristics for a 
constant spanwise pressure gradient corresponding to AC p = 0.15. Also 

c 

plotted in figure 19 for comparison are the results for N.F. and S.F. 
transition character!' sties on the laminar-flow glove for a spanwise pressure 
gradient of zero. It can be seen from these results that the Reynolds number 
at which S.F. triggers transition increases from 12 to 20 million as the 
spanwise pressure gradient is changed from zero to AC p = 0.15 (see figure 

18(a)). Furthermore, notice the change in N.F. transition characteristics for 
these two cases. Calculations indicate that turbulent separation does not 
occur in either of the cases shown in figure 19. 

Figure 20 shows the plots of C D versus freestream Reynolds number for 

wake 

the various constant spanwise gradients, including zero, considered in the 
present study. Thus, the results shown in figure 20 illustrate the ability of 
the present theory for increasing the design Reynolds number for minimum wake 
drag by controlling spanwise pressure gradients. It is necessary to point out 
here that constant spanwise pressure gradient increases the shock strength of 
the inboard spanwise location, n = 0.49 and decreases it at the outboard 
station, n = 0.91. It is likely that these shock strengths (figure 18(a)) 
will yield different pressure recoveries than that assumed. 

Another way to obtain a spanwise gradient is to fix the shock strength 
and vary the pressure levels near the leading edge. Figure 21 shows the 
computational results for transition due to N.F. and S.F. phenomena for such a 
pressure distribution i.e., the linearly varying spanwise pressure gradient, 
shown in figure 18(b) (aC r = 0.15) . Also shown in figure 21 are the 

theoretical transition characteristics for an X-29 glove when the spanwise 
pressure gradient is assumed to be zero (infinite-span swept wing 
conditions). Results of figure 21 illustrate that the triggering of 
transition due to S.F. is delayed from a Reynolds number of 12 million to 16 
million due to the effect of spanwise pressure gradient. In addition, the 
shock strength remains constant along the span when the linearly varying 
spanwise pressure gradient is assumed. Thus, when the linearly varying 
spanwise pressure gradient is specified, the shock wave drag remains the same 
as for a zero spanwise pressure gradient, and also the adverse effect of 
increased shock strength on the laminar glove design is eliminated. 
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versus freestream Reynolds 


Figure 22 shows the plots of computed C D 

wake 

number for linearly varying spanwise pressure gradients for several values 
of AC P . Also shown are computations for a spanwise pressure gradient of 
*1 

zero for comparison purposes. The results shown in figure 22 further 
illustrate how the present theory can be used for delaying transition Reynolds 
number due to triggering of S.F. instabilities by approximately 6 million 
above the zero spanwise-pressure gradient value by the use linearly varying 
spanwise pressure gradients. Neither induced drag, Cp., or shock wave drag 
are adversely affected. 1 


Comparison Between Equivalent Forward and Aft 
Swept Wing Concepts for Laminar Flow 

A comparison of planforms and local sweep variations for FSW and ASW with 
an aspect ratio of 4 and a taper ratio of 0.4 is shown in figure 11. The 
pressure distribution and spanwise pressure gradient selected for making a 
comparison of laminar flow regions on equivalent FSW and ASW is a linearly 
varying spanwise pressure gradients with aC p = 0.15. This spanwise pressure 

*1 

gradient corresponds approximately to that provied by an elliptical spanwise 
load distribution for both FSW and equivalent ASW. The results of 
computations for transition locations as a function of the freestream Reynolds 
number for these "equivalent" wings are shown in figure 23(a) and calculations 

of C n versus freestream Reynolds number are shown in figure 23(b). Figure 
wake 

23(a) shows that extent of laminar flow for FSW is 15 to 35 percent more than 
the equivalent ASW for Reynolds numbers from 8 to 30 million. The wake drag 
for FSW at freestream Reynolds number of 16 million, which corresponds to 
flight condition for the X-29 aircraft at = 0.905 and = 0.435, is lower 
by 21 counts than the equivalent ASW. 

The comparison shown in figures 23(c) and 23(d) assumes a favorable 
spanwise pressure gradient that provides a spanwise flow that is opposite to 
that normally obtained on ASW; this provides a spanload that deviates from the 
elliptical and causes an increase in C D . for ASW. Results of computations of 

the extent of laminar flow, shown in figure 23(c), indicate the extent of 
laminar flow for FSW is still larger by 7 to 20 percent over the "equivalent" 
ASW. A comparison of C n values under this assumption, as shown in figure 

23(d), indicates that for ASW the minimum value of C n occurs at a Reynolds 

u wake 

number of 11.5 million, whereas for FSW, C n becomes a minimum at a 

u wake 

freestream Reynolds number of 16.5 million. Moreover, values of C n for 

u wake 

the range of Reynolds number from 8 to 30 million are lower for FSW than the 
equivalent ASW. 
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It is necessary to point out that even though the FSW concept, in 
general, is better for a laminar-flow wing than the equivalent ASW, it is 
possible to determine other pressure distributions and spanwise pressure 
gradients for ASW that yield larger extents of laminar flow than shown in 
figure 23. The insight for obtaining the conditions for large laminar flow 
regions for highly swept ASW, especially for supersonic flow conditions, can 
be obtained by going through an optimizaiton process for chordwise pressure 
distribution and spanwise pressure gradients as has been done for FSW of X-29 
aircraft in the present studies. 


CONCLUDING REMARKS 

(1) The present integral methods for computing laminar boundary-layer 
properties and transition have been extended to take into account boundary- 
layer wash-in/wash-out phenomena and the effects of spanwise pressure gradient 
on transition for the purpose of delaying transition due to S.F. instabilities. 

(2) The above method has been applied to predict transition phenomena on 
the clean-up glove of the F-14 aircraft. Good agreement has been obtained 
between the predicted location of transition and experimental data as a 
function of Reynolds number, Mach number, and wing leading edge sweep. The 
results of the present studies for the F-14 aircraft wing glove, as well as 
those of references 5, 6, and 7, indicate that once transition occurs due to 
S.F., the transition location moves upstream very fast with even slight 
increase in freestream Reynolds number. Thus, by the use of the present 
method the wing-section geometry can be designed and spanwise pressure 
gradient specified in such a manner that transition due to S.F. pheonomena can 
be substantially delayed. This facilitates obtaining larger laminar-flow 
regions and making the wing performance less susceptible to minor changes in 
Reynolds number and Mach number. 

(3) Comparisons of the computed displacement thicknesses, momentum 
thicknesses, and form factors in turbulent boundary-layer regions on the glove 
of a F-14 aircraft by the present turbulent boundary-layer method with 
experimental data for several freestream conditions and wing sweeps indicated 
good agreement. This implies that not only is the present turbulent boundary- 
layer theory quite reliable, but the transition location prediction is 
accurate as wel 1 . 

(4) The extent of laminar flow computed for the F-14 wing glove for a 
leading-edge sweep of 25° and M^ = 0.7 was larger than that computed for the 
smaller sweep value of 20°. The experimental data obtained from hot film 
traces also revealed this unusual pheonomena. This provides additional 
evidence of the ability of the present method to optimize wing sections as 
well as the planforms for obtaining large extents of laminar flow. 
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(5) The effective sweep procedure which is normally used for ASW, when 
used to compute the streamwise pressure distribution on the FSW of the X-29 
aircraft, did not yield good correlation with experimental results. The new 
effective sweep procedure established during the present study for FSW gives 
reasonable predictions of streamwise pressure distributions when compared to 
experimental data for the X-29 wing. 

(6) By the use of the method of characteristics, in conjunction with the 
effective sweep procedure established for FSW, section geometries were 
established for the outboard glove for the wing of the X-29 so as to obtain 
the desired chordwise pressure distributions for leading-edge sweeps of 30° 
and 25°. The computed transition and wake drag characteristics were also 
computed and compared for these sweeps. The methods were also used to compute 
transition and wake drag for the basic outboard wing section. This was done 
for assumed "infinite-span swept-wing" conditions. The wake drag for the 
investigated glove section was quite superior to the basic section, and the 

glove with a leading-edge sweep of 25° had a minimum value of C n at 

wake 

freestream Reynolds number of 16 million as compared to 10.5 million for a 
sweep of 30°. 

(7) Several spanwise pressure gradients, both constant and linearly 
varying, for x/c = 0 to x/c = 0.6 on the upper surface gloves, were 
investigated for their ability to delay the transition due to S.F. phenomena 
at = 0.905. It was concluded from these parametric studies that both 
linearly varying and constant-increment spanwise pressure gradients provided 
beneficial effects on skin friction drag at the flight condition considered. 


(8) Computational experiments were performed to determine the extent of 
laminar flow on FSW and "equivalent" ASW for identical chordwise pressure 
distributions and spanwise pressure gradients. This comparison indicated that 
the FSW concept is superior for achieving large regions of laminar flow and 
lower wake drag than an "equivalent" ASW concept. 
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APPENDIX A 




Figure Al. - Coordinate system used for wake drag derivations. 


Consider control volume ABCD-EFGH as shown in figure Al above for a wing 
with sweep angle A. The x-coordinate is perpendicular to the leading edge, 
the z-coordinate is parallel to leading edge, the y-coordinate is 
perpendicular to the wing surface, S is equal to wing span as shown, and c 
is the chord of wing in the direction of flow. 

Net mass balance °° 

across the control = / p^U^S cos A dy - / puS dy = 0 (Al) 

volume -« -» 


Drag = Net momentum flux across control volume 


2 

= / plr S cos A dy - / puS(u cos A + w sin A)dy (A2) 


Multiply equation (Al) by U m and subtract from (A2) 
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2 2 

Drag = / pU uS(cos A + sin A)dy - / puS(u cos A + w sin A)dy 


= S cos A p / pU (1 - — )dy 

K ooooJ pu v u } J 


— OO 00 00 


+ S sin a p u w 7 (1 - £-)dy 

CO 00 00 ^ pU v W y 


—00 00 OO 


or 


. Drag 

wake 75 - p U S*c cos A 

C oo oo 

Drag 

12 2 

- 5 - p u sec A*S*c cos A 
^ 00 00 


where 


0 


uu 


0 , 


uw 


— r OO / 

2{ — — cos A + — 
1 c c 


2 

sin a} 


O = f (1 _ M-)dy 

UU j p u ^ u 


00 — 00 00 00 


0 = / _pJL_ (1 . w_) dy 

UW J p U ^ W ' J 

00 —00 00 00 00 


(A3) 


(A4) 


For swept wings as well as bodies of revolution the boundary-layer 
parameters at the trailing edge can be calculated by the theoretical methods 
described in the main body of this paper. Equation (A4) gives C D in 
terms of conditions downstream at infinity. The following derivatY&kf show 

how boundary-layer parameters at the trailing edge are related to viscous flow 
quantities at infinity. 

The momentum integral equation in the wake in the streamwise direction, 
either for swept wing or for body of revolution can be written as 


where 




du 

dx 


6 *x 1 dp * 

(2 + — ) + — 1 — 

1 e x °e dx 


e x = ° 


©X = 0 UU of equation (A4) 


(A5) 
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By integrating equation (A5) in the wake from the trailing edge to downstream 
infinity, and after algebraic manipulations the following equation is obtained 


[*"( e x P e (F“) 


(H„+2) - » u 

)] ■ / 

T.E. T.E . 


(A6) 


The hypothesis is made that in the wake the following relationship is valid. 


tn 



tn 


'T.E. 



H 





(A7) 


The following equation can be derived by making use of Stewartson's trans- 
formation for a Prandtl number equal to one. 


H = H t 
x T 


(i * ^ 


(A8) 


H x = streamwise form factor in physical plane 

H = streamwise form factor in transformed plane 
X T 

Substitution of equations (A7) and (A8) into (A5) and simplification results 
in the following equation 


uu 


(») 


= 0 


uu 


T.E. 


(1 


^M 2 


e,s 


..-2.5 

■’T.E. 


(1 


r-1 


M 2 )- 2 * 5 

oo ' 



(A9) 


where. 


P = {H 


T.E. 


+ 5 + 


( y - Dm !}/2 


m Root mean square Mach number in the wake from trailing 
e edge to infinity downstream 


Equation (A4) for C n for a swept wing contains the crossflow momentum 

u wake 

thickness term, 0 . In order to evaluate this crossflow component of 

uw 

oo 

momentum thickness, use is made of the momentum equation in the spanwise 
direction for an infinite-span swept wing. Thus, 
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(A10) 


pu 


aw 


+ pv |^ + p w |^ = - IE + JL r„ 

H 3V M 37 37 3v 9 yJ 


ax FT ay ' az az ' ay 
The infinite-span swept-wing assumption implies 


h < 1 -° 


and the continuity equation under the same assumption yields. 


(All) 


p v=-/ (pu)dy 


0 


(A12 ) 


Along the locus of minimum velocity in the wake |y = 0. Substituting (All) and 

(A12) into (A10) and integrating with respect to y in the wake from 0 to 
6 Z , the following form of momentum integral equation in the wake can be 
deri ved. 




e e 


37 [ p e u e w e e J = 0 


0 = ( 
UW *■ 


p e u e w fc 
T.E.w e T.E.w e T . E . 


u ^ w ^ 0 uw^T.E. 


(A13) 


Thus, C D can be written from equations (A4), (A9), and (A13) as: 
wake 


c [ (0 uu } T.E. (p^T.E. (ir) cos2a T.E. 


wake 


w 


+ ( 0 uwH.E. (p^T.E. T " E ') sln ^ A T.E . J 


(A14 ) 


Equation (A14) contain terms evaluated at the trailing edge, such as 
^uw^T.E.’ ^ 0 uuH.E.’ ^ p e^T.E.’ an< * ^ u e^T.E ’ These values are computed in 
the framework of computer programs for boundary-layer methods, (refs. 5, 6, 7, 
and 8) and pressure-distribution prediction methods such as those of references 
9 and 10. 
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M co sin A 



u(s,^ )— Velocity profile in 

streamwise direction 

v(s, ^ ) Velocity profile in 

direction normal to wall 

w(s,^ )— Velocity profile in 
spanwise direction 


Figure 1. - Definition of coordinate axis system. 


P + 



dP 





Figure 2. - Schematic of control volume used to account for boundary-layer 
wash-in/wash-out phenomena. 



a o e P,n /dM e \ 1 

av m \ dTl 7(1 + -2M 2)2 

Figure 3. - Criteria for determining laminar boundary- layer neutral 
instability. 



Figure 4. - Criteria for determining laminar boundary-layer transition 




Free-stream turbulence, % 

Figure 5. - Effect of free-stream turbulence on the difference of transition 
and instability Reynolds numbers. 



• Three rows of flush static pressure orifices 

• Three rows of surface pitot tubes 

• Five hot-film gauges 

• Two boundary-layer rakes 


Figure 6. - Planform of an instrumentation on the gloved F - 1 4 wing. 
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(Moq^ 0.700, A L e = 20.0°, Rco = 13 x 1q6 ) 


11 = 0.52 


il = 0.67 


11 = 0.84 




1.0 


(a) Experimental pressure distributions and hot-film traces. 


(M jo = 0.697, A LE = 20°) 



(b) Computed transition locations as a function of free-stream chord Reynolds 
number. 

Figure 7. - Results for the F-14 "clean-up" gloved wing for a leading- edge 
sweep of 20° and of 0.7. 
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(Mqj = 0 . 697 , A le= 20 ° , Roo = 13.6 x 10 6 ) 



(c) Comparison between computed and flight boundary- layer parameters for 
= 13 x 10°. 

Figure 7. - Concluded. 




I 

E 
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Tl = 0.52 


ri = 0.67 


0 - 0.84 




HOT FILM TRACES 

Station 3 


(a) Experimental pressure distributions and hot-film traces for 
R « 11 x 10 6 . 

GO 




HOT FILM TRACES 

Station 2 


(b) Experimental pressure distributions and hot-film traces for 
R m = 23 x 10 6 . 

oo 

Figure 8. - Results for the F- 14 "clean-up" gloved wing for a leading-edge 
sweep of 20° and of 0.8. 
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NF (Normal flow) 


Transition 


SF (Spanwise flow) 


Flight 8 , Station 3 


Flight 18, Station 2 


Instabilit y 


Roo/10 6 

(c) Computed transition and instability locations versus free-stream Reynolds 
number. 


C D wake 


(A le = 20° > Station 3) 


(d) Computed C n versus free-stream chord Reynolds number 
u wake 

Figure 8. - Continued. 





0 .2 A .6 .8 1.0 


X/C 

(e) Comparison between computed and flight boundary-layer parameters for 
R ro = 12.4 x 10 6 . 

Figure 8. - Concluded. 
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HOT FILM TRACES 

Station 3 


(a) Experimental pressure distributions and hot-film traces. 



(b) Computed transition locations as a function of free-stream chord Reynolds 
number*- - 

Figure 9. - Results for the F- 14 "clean-up" gloved wing for a leading-edge 
sweep of 25° and of 0.7. 



0/c, 6 



(c) Comparison between computed and flight boundary-layer parameters for 

R = 10.8 x 10 6 . 

00 


Figure 9. - Concluded. 


T1 = 0.52 


n = 0.67 


il = 0.84 



C 0.1 0.2 0.3 0.4 0.5 

HOT FILM TRACES 


(a) Experimental pressure distributions and hot-film traces for 
R = 9.5 x 10 6 . 

oo 


t ) = 0.52 r| = 0.67 n = 0.84 



c 0.1 0.2 0.3 0.4 0.5 

HOT FILM TRACES 


(b) Experimental pressure distributions and hot-film traces for 
R^ = 19.2 x 10 6 . 

Figure 10. - Results for the F-14 "clean-up" gloved wing for a leading-edge 
sweep of 30° and M of 0.7. 




R co /io® 


(c) Computed transition and instability locations versus free-stream Reynolds 
number. 

Figure 10. - Concluded. 
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(b) Local sweep variations with respect to x/c. 

Figure 11. - Comparison between "equivalent" forward and aft swept wings. 
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T1 



Figure 12. - Planform of X-29 aircraft. 



(Z/C) fg = (Z/C) 5 /cos A upper and lower for Aft Sweep (N = normal, S = streamwise) 

(Z/C) n, upper = (^) s, upper ^ Forward Sweep 

(Z/C) N) , ower = (Z/C) s, upper * ( 1/cosA ^ + < Z/C) S, lower * ( 1/cos A * 

(M L ) N = (M l ) s • cos A x /C =0.25 

(C p ) streamwise = < C P> N ‘ cos 2 A where A = A LE + < A TE ‘ A LE ) * ( x/c > 

(Cl) 2-D = (C L ) s /cos 2 A X/C=0.25 

Figure 13. - Sweep procedure for forward- and aft-swept wings. 


} Forward and 
Aft Sweep 
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r r ^° rwar d-swept 

■■ 6 Jf x °i JD_ — O— 

J \ 


Aft swept 

- Experimental upper 


C p -.2 


□ □ □ □ 


Experimental lower 


0 1 -2 .3 .4 .5 .6 .7 .8 .9 1.0 

X/C 

Figure 14. - Computational results for the pressure distributions on the 
X-29 wing section situated on both forward and aft winqs 
/u m = 0.91, r, = 0.70, a le = 30°' ’ 



C D wake 
10-4 



(b) C n for baseline versus R and comparison with NLF-Glove. 
u wake 

Fiqure 15. - Computed results for X - 29 baseline section for rj = 0.7, 

Moo - 0-91- 
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(c) Transition characteristics for leading-edge sweeps of 30° and 25 



Figure 17. - Concluded. 



C P, upper 


C P, upper 


Figure 18. - 




(b) Linearly varying spanwise pressure gradient. 

Schematics showing types of spanwise pressure gradients studied. 
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X/C .3 


(M co =0.905, A le =30°) 

NF,ACp = .15 
r c 

NF (Normal flow), ACp = 0 ^ 

c \ 


SF, ACp = .15 —A 
C \ 


SF (Spanwise flow), ACp =0 

.1 — c 


Re ©/™ 6 


Figure 19. - Comparison of the effect of a constant spanwise pressure 
gradient of aC p = 0.15 on transition with that for zero 
,. c 

spanwise gradient, i.e., with aC p = 0. 

r c 


(M „,= 0.905, A LE =30°) 


AC p = 0 




versus with ACp as parameter for X- 29 glove 












NF (Normal flow), ASW 


NF,FSW 


SF (Spanwise flow), FSW 


SF, ASW 


ASW - Aft Swept Wing A j_j= = 45° 

■1 “ FSW - Forward Swept Wing A j_£ = 30° 


R^/IO 6 

(c) Transition location versus for ACp = -0.15. 


Forward Swept Wing A le ^ 30° 
44 f Aft Swept Wing A LE = 45 ° 


versus R^, for aC p 
ke r l 

Figure 23. - Concluded 



